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Abstract. We consider a class of unstable surface growth models, d t z — —d x J, developing a mound 
structure of size A and displaying a perpetual coarsening process, i.e. an endless increase in time of A. The 
coarsening exponents n, defined by the growth law of the mound size A with time, A ~ t n , were previously 
found by numerical integration of the growth equations [A. Torcini and P. Politi, Eur. Phys. J. B 25, 519 
(2002)]. Recent analytical work now allows to interpret such findings as finite time effective exponents. The 
asymptotic exponents are shown to appear at so large time that cannot be reached by direct integration 
of the growth equations. The reason for the appearance of effective exponents is clearly identified. 

PACS. 68 Surfaces and interfaces - 81.10.Aj Theory and models of crystal growth - 02.30.Jr Partial 
differential equations 

1 Introduction in the surface current J Possible noise sources, shot 

noise first, will be neglected. 



H . 

' I n this manuscript we are interested in studying a class of 
one dimensional growth equations, having the conserved 

form fa = F - d x J : the dynamics of the local height The current j may have a plethora of different forms, 
z{x,t) of the surface is determined, apart from a trivial depe nding on the details of the atomistic processes 0. 
constant term F describing the deposition flux, by the For Molccular Beam Epitaxy gl|S], a widely used tech- 
processes occuring at the surface, which are all included nique for growing metal and semiconductor thin films with 
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the form The study of the stationary periodic solutions is of 

great importance for the dynamics as well. In Ref. 0, 



d t z = -d x (u xx + j(u)) (1) 

where the constant term F has been included in the left 
hand side by redefining z — z — Ft, and u = z x = d x z is 
the slope of the surface. It is worth noting that taking the 
spatial derivative of both sides, we get Ut = ~d xx {u xx + 
j(u)), i.e. a generalized Cahn-Hilliard equation jHj. With- 
out the double derivative (— d xx ) we get the corresponding 



for Eqs. IL'l-il) and other classes of models, it has been 
shown that the surface undergoes a coarsening process if 
and only if dX/dA > 0. In simple words, the wavelength 
of the mound structure (emerging from the linear insta- 
bility) increases in time if the wavelength A (A) of the pe- 
riodic steady state increases with the amplitude A. Even 

more importantly, the knowledge of the stationary peri- 
non-conserved models, called generalized real Ginzburg- ,. . ,. , , _ ;1 , 

odic solutions allows to determine 6 the coarsening law 

Landau (or Allen-Cahn) equation 6 . In summary, we are , „ , ■ r i \ 

X{t) (see the next Section for more details). 

considering the two classes of equations 

u t = -d xx (u xx +j(u)) conserved (2) 2 Numerics vs Analytics 

u t = u xx + j(u) non-conserved (3) Jn thc following we are focusing on a class of models de- 

The linear stability analysis of the flat interface, z = fined bv the currents 

zq + eexp(u>t + iqx), can be easily applied to Eqs. j2EB j(u) = - , (6) 

(1 + u 2 ) a 

gi ym g 

which correspond to the potentials 

w = j'{0)q 2 - q A conserved (4) ( , _ __L 1 m 

K ) 2(a-l) (l + u 2 )"" 1 ' K ) 

lu = j'(0) — q 2 non-conserved, (5) 

They were introduced in a previous paper on this jour- 

so that an instability appears (i.e. u)(q) > for some q) nal 8 and were called a— models. 

if j'(u = 0) > 0. Without loss of generality, we may as- It is straightforward to check that V(u) has a mini- 

sume j (0) = 1. The steady states are determined by the mum in u — and goes to zero for large u, as V(u) ~ 

equation u xx = —j(u) both in the conserved (J2J) and non- — l/|w| 2 ( a_1 '. The curves X(A) can be found numerically, 

conserved © case. The solutions u(x) correspond to the but the sign of dX/dA can be easily deduced from the be- 

trajectories of a fictitious particle moving in the symmet- havior of V(u) at small and large u. At small u, V(u) « 

ric potential V(u) — J duj(u) — \u 2 + higher order terms. V(0) + ^u 2 — |a 4 , so that the quartic term is negative. At 

The oscillations in the potential well correspond to peri- large u, V(u) increases and goes to a constant value. Both 

odic steady states of the variable u, with wavelength A these features are signatures [I0J for a positive dX/dA, i.e. 

and amplitude A. for a perpetual coarsening. 
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2.1 Old results diffusion coefficient, which describes the dynamics of the 

local phase, when the periodic stationary solution is per- 

In Ref. |Hj we studied the coarsening exponent n (A(t) ~ 

turbed. This approach is applicable to large classes of 

t n ) for the a— models. The non-conserved version, Eqs. ilol(il) . 

models, both conserved and non-conserved. As for the 

allowed for an analytical treatment which followed an ap- 
proach due to Langer f5j: it consists in evaluating the 



n 



—models, X(t) is deduced from the relations [Hj 



most unstable eigenvalue of the linear operator describ- iOjj ^ 1 conserved models (TO) 

IX' {A) ~ t conserve m0 es I J 
ing perturbations of the periodic stationary solution. This j(A) 1 

- - non-conserved models (11) 



method gave the results 



J A' ( A) t 



2 



I a < 2 In Eqs. (|1C)I1 1|> . all the quantities refer to the periodic 

non-conserved (Ref. j^j) (8) 
a > 2 steady states, u(x + A) = u(x), satisfying the equation 



3a-2 

The conserved version, Eqs. I|2I6|) . of a— models did not u xx + j(u) = 0: A is the amplitude (the maximal pos- 

allow for an equally rigorous approach. Therefore, we inte- itive value of u(x)); X is the oscillation period (i.e. the 

greted numerically [8] the growth equations d t z = —d x J: wavelength) and X'(A) is its derivative with respect to A; 

it appeared that coarsening exponents agreed fairly well J is the action variable, defined by J = § dxu x ; finally, 

with the relations / = § dxu 2 . 

j a < 2 For large A, we can split the motion of the fictitious 
conserved (Ref. 0) (9) 

a > 2 particle in the potential V(u) in a region close to the ori- 

These results have been interpreted by doing the fol- gill; | u | < ^ and in the complementary regions A < 

lowing ansatz: "as for the coarsening exponent, passing | u | < A Aq is choosen so that in the reg ions |u| > A , j(u) 

from the non-conserved to the conserved models is equwa- and V (u) can be approximated by their asymptotic ex- 

lent to replace (-d xx ) in Eq. 0) with 1/A 2 ." This recipe pre ssions j(u) ~ 1/lul 2 "" 1 and V(u) ~ - 2[hT)\ u \~ 2{a ~ 1] ■ 

allows to get the conserved coarsening exponents © from gince y{A) goeg tQ & constant for diver g ing A? thc 

the non-conserved ones © in a straightforward manner. motion of ^ partidc m thfi gmaU u region < Aq) 

The previous picture appeared to be reasonably correct does not depend Qn A There fore, every integral quan- 

until a more general analytical approach has been re- tity ^ j J} ig thfi gum of a constant ternl; coming from 

cently developed. integration in the small u region, and an asymptotic 



2.2 New results 



A— dependent term, coming from the integration in the 
large u region. In all cases (with one exception, see below) 
This new theory relies on the observation that the coars- the asymptotic contribution diverges with A and therefore 
ening law X{t) can be extracted from the so-called phase dominates. Such contribution can be evaluated by dimen- 
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sional arguments and, more rigorously, using the law of a 
mechanical similarity For example, X(A) can be sim- r ^,j 
ply deduced equating the 'acceleration' A/X 2 to the 'force' 
j(A) ~ l/A 2 " -1 , therefore getting A ~ A a . Similarly, we 
get / ~ XA 2 . 

As for J, the asymptotic contribution amounts to X(A/XY 
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A 2 - a , which diverges for a < 2 only. For a > 2 (this 10 10 10 10 10 10 
is the exception mentioned above) the asymptotic contri- 
bution of J vanishes, indicating that the small u region, 
giving a constant contribution, dominates. In conclusion, 
J ~ A 2 - a for a < 2 and J ~ 1 for a > 2. 

If we replace the previous relations in Eqs. I|1UI11|) . we 
get Eq. (JHJ) for the non-conserved models, but we get a 
constant coarsening exponent, n = \, for the conserved 
models, in sharp contrast with Eq. (0. In conclusion, 
the Langer-type approach and the recent theory based on 
the phase diffusion coefficient give the same results © 
for the non-conserved a— models. For the conserved mod- 
els, where the Langer-type approach is not applicable, the 
recent theory seems not to agree with the numerical re- 
sults (0 found in Ref. 8; via direct numerical integration. 
This disagreement also calls for reconsidering the ansatz 
(—d X x) 1/A 2 . Next Section is devoted to understand 
the origins of this discrepancy (therefore, we will limit to 
a > 2). 
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10 
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Fig. 1. The coarsening law X(t) for the conserved model a = 3, 
as derived by Eq. 11100 after numerical determination of the 
steady states u(x). Full and dashed lines refer to the asymptotic 
fit and to the fit in the region A « 10 2 , respectively. Inset: The 
quantity I/X 2 as a function of A, for the same conserved model. 

performed in Ref. [E] up to A « 10 2 . A numerical fit in this 
region (dashed line) gives an effective esponent n — 0.231, 
which is in perfect agreement with such simulations and 
with the relation n = a/ (5a — 2). Conversely, a fit in the 
asymptotic region (full line) gives n — 0.25, as expected. 
The numerical results found in [S] are therefore interpreted 
as finite-time exponents: formulas @ are applicable up to 
A w 10 2 , but the asymptotic exponent for the conserved 
model is j. 

Let us now discuss the origin of such finite-time ex- 
ponents. As discussed in the previous Section, formulas 
(0 are correct insofar as the ansatz (— d xx ) — * 1/A 2 is 
correct. Its validity corresponds to say that 



3 The origin of the effective exponents (7) , ~ ti ( 7 ) • ( 12 ) 

\ 1 / conserved A \ 1 / non-conserved 
First of all, let us determine numerically X{t) from Eq. JTUJ If we focus on our exact relations UiiHU. this rclation 
at all times. Results for a — 3 are shown as circles in Fig. 1 would imply 

(main). The direct integration of the growth equation was /(S^A) A 2 J(8aX) ' 
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current 10 and termed a— models, all display perpetual 
coarsening, X(t) ~ t n . 

For the non-conserved models, n depends on a accord- 
ing to formula ijSJ): therefore, the theory based on the phase 
diffusion coefficient |Hj confirm previous results [8] based 
on a Langer-type approach [§]. For the conserved models, 
recent theoretical results [H] give a constant coarsening 
exponent, n — j, at odd with our previous numerical re- 
sults, Eq. ©. We have explained the effective exponent 
n = a/(ha — 2) as a finite-time exponent, whose appear- 
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Fig. 2. The scaling of the maximal force j(A) (left) and of the 
quantity A'(^4) (right), as functions of A. 



i.e. I/X 2 ~ J. Since J is constant for a > 2, the formula ance is due to thc f act that tne quant it y //A 2 , instead 

n = a/(5a - 2) would be correct if (I/X 2 ) would be con- f decreasing as a power law, is approximately constant 

stant as well. In the inset of Fig. 1 we plot the numerical f or not too large A (see Fig. I, inset). This constant be- 

results for I/X 2 vs A: for very large A, I/X 2 ~ A" 1 / 3 de- haviour is equivalent to assume that the operator (-d xx ) 



creases, but for A ~ I0 2 it has a maximum. In other words 
in the region of wavelengths which can be reasonably in- evaluated 
vestigated with the direct integration of the growth equa- 
tion, the quantity I/X 2 is approximately constant, which 
implies n a/ (5a — 2). 

In order to support the idea that the origin of the effec- 



in Eq. J5J can be effectively replaced by I/A 2 , when n is 



It is worthnoting that the result n = j for the con- 
served models was firstly found by Golubovic using 
some dimensional arguments (see Section 7 of Ref. jS] for 
more details). The main problem with this approach is 



tive exponents is indeed the maximum in I/X 2 , in Fig. 2 that y. algo giyes & constant coarsening expone nt n = \ for 
we also plot the two other relevant quantities appearing in ^ non _ conscrvcd modelS) which is wrong , Thc rcason of 
Eq. HDD, J (A) and dX/dA: none of them has any special tUg failure ig ^ from SectionE3 because dimensional 
behaviour for small A. analysis can be valid only if the A-dependent asymptotic 

contribution to the integral quantities A, /, J outnumbers 
. _ the contribution from the small u region. This is not thc 

4 Conclusions 

case for the non-conserved models and a > 2: exactly the 
In this short note we have reconsidered a class of conserved class of equations where dimensional analysis fails. 
J21 and non-conserved © growth models, in the light of Similar arguments can be used to understand the fail- 



recent theoretical results (Hj- These models, defined by the ure of the ansatz (— d xx — » 1/A 2 ), which has clearly to 
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do with dimensional analysis. The weak point is that this time step up to four times, while the integration scheme 

ansatz is applied along with the correct results for the becomes rapidly unstable by considering a coarser space 

non-conserved models, where dimensional analysis works grid. As a matter of fact, we cannot expect to lower the 

for a < 2 only: therefore, the ansatz gives the final correct CPU time more than a factor ~ 10, which renders still un- 

result (for the conserved models) for a < 2 only. feasible the observation of the asymptotic exponents (by 

It is also worth stressing that a more rigorous applica- the wa ^ even a lowering up to a factor 10 5 would make it 

tion of dimensional arguments 13 to the conserved mod- unfeasible), 
els gives an inequality, n < 4, which does not allow to 

discriminate between n = \ and n = a /(5a -2). References 
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